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Abstract 

The R-Parity symmetry Violating (RPV) version of the Next-to-Minimal Supersymmetric Stan- 
dard Model (NMSSM) is attractive simultaneously with regard to the so-called /i-problem and the 
accommodation of three-flavor neutrino data at tree level. In this context, we show here that if 
the Lightest Supersymmetric Particle (LSP) is the gravitino, it possesses a lifetime larger than the 
age of the universe since its RPV induced decay channels are suppressed by the weak gravitational 
strength. This conclusion holds if one considers gravitino masses ~ 10^ GeV like in supergravity 
scenarios, and is robust if the lightest pseudoscalar Higgs field is as light as ~ 10 GeV [as may 
occur in the NMSSM]. For these models predicting in particular an RPV neutrino-photino mixing, 
the gravitino lifetime exceeds the age of the universe by two orders of magnitude. However, we find 
that the gravitino cannot constitute a viable dark matter candidate since its too large RPV decay 
widths would then conflict with the flux data of last indirect detection experiments. The cases of a 
sneutrino LSP or a neutralino LSP as well as the more promising gauge-mediated supersymmetry 
breaking scenario are also discussed. Both the one-flavor simplification hypothesis and the realistic 
scenario of three neutrino flavors are analyzed. We have modified the NMHDECAY program to 
extend the neutralino mass matrix to the present framework. 
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I Introduction 



In supersymmetry, the superpartner of the graviton, namely the so-called gravitino, plays a central 
theoretical role as it constitutes the gauge fermion of supergravity theories [H [2] . From the cosmological 
point of view, the gravitino may be the Lightest Supersymmetric Particle (LSP), depending on the 
way supersymmetry is broken, and hence constitute a stable dark matter candidate p]. Interestingly, 
one could even assume supersymmetric models where the R-parity symmetry [4] is violated. Then the 
gravitino becomes unstable due to new decay channels into Standard Model (SM) particles induced by 
the R-Parity Violating (RPV) interactions. However, because of the small gravitational interaction 
and possibly weak RPV couplings, the gravitino lifetime can still exceed the age of the universe (by 
several orders of magnitude) in the particular case of RPV bilinear terms [5] remaining thus a realistic 
dark matter candidate. This scenario can yield a complete cosmological framework also incorporating 
leptogenesis and nucleosynthesis, as recently shown in Ref. [6j. In such a scenario, the gravitino decays 
into SM particles lead to specific signatures in high-energy cosmic rays; the produced flux of gamma 
rays and positrons can even account [7] respectively for the extragalactic component of the excess in 
the HEAT [8j data. The decaying gravitinos could also constitute an interpretation (see below) of the 
exotic positron source recently discovered by the PAMELA Collaboration [9]. 

Another advantage of the existence of RPV interactions is the induced mixing between left-handed 
neutrinos and neutral gauginos. This mixing mechanism generates neutrino masses economically, 
without extending the SM field content to additional right-handed neutrino fields. There exists a 
tension between generating a sufficiently large neutrino Majorana mass scale [requiring strong RPV] 
and, at the same time, keeping the gravitino lifetime larger than the age of the universe [imposing 
weak RPV]. Nevertheless, this tension leaves some acceptable windows in the parameter space [5] 
for example with gravitino masses around 100 GeV in the Minimal Supersymmetric Standard Model 
(MSSM). 

In the present paper, we re-consider the above scenario in the framework of the Next-to-Minimal 
Supersymmetric Standard Model (NMSSM) (see recent reviews |11] for the phenomenological 
studies). An increasingly important virtue of the NMSSM is that it improves the 'little hierarchy' 
problem originating from the requirement of large soft supersymmetry breaking masses compared to 
the ElectroWeak (EW) scale. It also provides a solution to the so-called /i-problem by arranging the 
vacuum expectation value (vev) of a new gauge singlet scalar field of order of the supersymmetry 
breaking scale, so that the n parameter turns out to be at the EW scale. We have shown in a previous 
work [12] that in the NMSSM with bilinear RPV terms, two non-degenerate massive neutrino states 
can emerge at tree level, in contrast with the MSSM case where only one neutrino eigenstate acquires 
a mass (at tree level). Hence, the global three-flavor neutrino data can be accommodated [13j, at tree 
level, without appealing to loop dynamics which is vulnerable to model-dependent uncertainties. 

In such a framework, both the gaugino mass matrix and RPV mixing terms are modified w.r.t. the 
MSSM case. Furthermore, the NMSSM provides the possibility0of reducing significantly w.r.t. MSSM 
the mass of the lightest pseudoscalar Higgs field [which can constitute partially the gravitino decay 
final state]. Moreover, the specific gauge singlet scalar field gives rise to new contributions to the 
gravitino RPV decay. Hence, the double question in the NMSSM on the possibility of generating 

^Such a possibility was in particular motivated by interpreting the well-known LEP excess at 2.3a in the Z + 2b 
channel - via the production of a light Higgs boson {nih — 99 GeV) [14] - but this NMSSM interpretation seems to be 
excluded by a recent ALEPH analysis [15| . 
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the correct neutrino mass scale, while still keeping the gravitino as a good dark matter candidate, is 
relevant and well motivated. 

We find in this paper that the answer concerning the gravitino stability can be positive for parameters 
passing the theoretical and phenomenological constraints implemented in the NMHDECAY code [16] 
and in particular for gravitino masses as large as m3/2 = 0(10^) GeV [tending to increase the total 
gravitino width]. This is satisfactory as it corresponds to the typical scale of supersymmetry breaking 
in supergravity scenarios. In both scenarios where the supersymmetry breaking is mediated purely by 
gravity [H [2] or partially by gravity and gauge interactions [T7] (such classes of hybrid models have 
been recently motivated in string inspired constructions [l^), the gravitino is the LSP in wide regions 
of the parameter space and its typical mass is in the range between 100 GeV and 1 TeV. Besides, 
thermal leptogenesis and universal boundary conditions for gaugino masses at the Grand Unification 
Theory (GUT) scale restrict m-^i2 < 600 GeV [19] and, on the other hand, the interpretation of the 
PAMELA positron anomaly imposes ^ 200 GeV |20j . 

Nevertheless, our second result on a gravitino LSP around ~ 10^ GeV is that its RPV decay rates 
lead to fluxes exceeding the total flux measured in last indirect detection experiments. 
We will finally discuss the pure case, still within the NMSSM, of Gauge-Mediated Supersymmetry 
Breaking (GMSB) where the LSP gravitino mass verifies typically ^ 1 sV [21j and the neutralino 
mass matrix receives some modifications [22\ . Conclusions here are more optimistic. 

Let us make some comments on the results obtained in the present paper. First, we have considered 
all the possible gravitino decays, namely into photon, Z, W bosons and (pseudo)scalar Higgs fields, 
in contrast with Ref. [5] where only the photon channel was considered. The new kinematically 
allowed channels open up due to the higher gravitino masses considered here. Secondly, the trilinear 
RPV couplings are also discussed. Furthermore, we investigate for the first time the realistic case of 
three flavors of neutrinos. We conclude that the complete RPV mixing obtained in this case does 
not invalidate the positive results obtained for the simultaneous solution of the neutrino mass and 
dark matter problems. Finally, the pure neutrino results include the various NMSSM-constraints 
implemented in NMHDECAY jl6j (that we have modified to include the neutrino components in 
the neutralino mass matrix) and are derived from a numerical treatment of the full neutral gaugino 
mass matrix, in contrast with preliminary work in Ref. |12j . The modified version of NMHDECAY 
thus contains an implementation of the whole neutralino mass matrix which includes, in the present 
context, some matrix elements induced by the RPV couplings and responsible for the mixing between 
the higgsinos and neutrinos (see Section . 

In addition, we will explore the alternative possibilities of the neutralino and sneutrino as stable LSP 
dark matter candidates, under the same assumption of the desired neutrino mass spectrum generated 
through the RPV mixing in the NMSSM. The philosophy is to establish a systematic list of the viable 
super symmetric dark matter candidates. Since this topic relies on the particle spectrum, we base our 
study on a systematic exploration of the parameter space. Our conclusions are also based on new 
calculations of the neutralino and sneutrino decay channels. 

Let us mention previous related works. The neutrino flux from direct gravitino decays has been 
computed in a simple scenario with bilinear R-parity breaking [23j . The diffuse gamma ray flux was also 
studied in Ref. [24j. A realization in minimal super symmetric left-right models within supergravity, 
with a gravitino LSP dark matter and R-parity breaking, was proposed in Ref. [25] (see also Ref. j26j). 
Finally, a scenario with R-parity violation in the right-handed neutrino sector was motivated by the 
PAMELA data (27j. 

Concerning supersymmetric models for neutrinos - independently of the dark matter problem - RPV 
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versions of the NMSSM have been previously studied in Ref. [28j, Ref. [29j and Ref. [30j . Besides, 
alternative super symmetric extensions with broken i?-symmetries have been proposed in order to 
address simultaneously the /i-problem and the neutrino mass aspect [311 [32l [33] (and see Ref. [M] for 
a gravitino dark matter decay study in this context). 

RPV supersymmetric scenarios have also been studied within the context of neutrino astrophysics (see 
e.g. [35i [36j ). Finally, for studies of gravitino dark matter without RPV couplings, see for instance 
Ref. [37| for the constrained MSSM or GMSB scenarios. 

The paper is organized as follows. In Section II, we elaborate the RPV scenario and describe the 
corresponding neutralino mass matrix. In the following sections, we discuss subsequently the cases of 
a sneutrino LSP (Section III), a neutralino LSP (IV) and a gravitino LSP (V). Finally, we conclude 
in Section VI. 

II The RPV version of the NMSSM 
II. 1 Theoretical framework 

We consider the NMSSM which possesses a superpotential containing two dimensionless couplings A 
and K in addition to the usual Yukawa couplings: 

Wnmssm = Y^^QiHJJ'r + YfjQiHdD^ + Vf^LiHaE^ + XSHuHa + ^kS^ (1) 

where Yj"''^'^ are the Yukawa coupling constants k are family indexes), and Qi, Li, Uf, Df, Ef, Hu, 
Hdi S respectively are the superfields for the quark doublets, lepton doublets, up-type anti-quarks, 
down-type anti-quarks, anti-leptons, up Higgs, down Higgs, extra singlet under the standard model 
gauge group. An effective /i term, given by X{s)HuHd, is generated via the vev of the scalar component 
s {{s)) of the singlet superfield S. 

In case of GMSB [22], there are additional terms that must be added to the NMSSM superpotential 
([!]). Those terms are given by: 

t^GMSB = ^FS + ^l'S^ (2) 

where ^i? (^u') is a new dimension-two (-one) parameter. These parameters in GMSB models are 
generated at low energy as S is coupled to the messenger sector. 

Recalling that there exists no deep theoretical principle in supersymmetry for the existence of an 
exact i?-parity symmetry [38], we adopt a generic approach by introducing both the bilinear and 
trilinear RPV terms characteristic of the NMSSM [the usual trilinear RPV couplings of the MSSM 
are considered in next section]: 

W = VFnmssm + fJ-iLiHu + XiSLiHu, (3) 
where /ij (Aj) are the dimension-one (dimensionless) RPV parameters. 

Within the usual NMSSM, only trilinear couplings with dimensionless parameters (like A and k) are 
kept in the superpotential, while dimensional parameters (like //) are generated from the vev (s). Here 
the RPV NMSSM superpotential ([3]), containing a ^iLiHu but no fiHuHd term, is assumed to arise 
in one of the scenarios proposed in [12j. 
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II. 2 Discussion on the trilinear RPV terms 



The most general RPV NMSSM superpotential also includes the renormalizable trilinear RPV inter- 
actions: 



It is remarkable that if the order of magnitude for Xij^ and X'-/^ is comparable to the values 
of the Yukawa coupling constants for the electron and the down quark [y^, Y"^], the exchange 
of leptons/sleptons and respectively quarks/squarks in one-loop processes can generate Majorana 
neutrino masses [39] in agreement with (or few orders below) the oscillation experiment results: 
10~^eV < mjy- < leV. This occurs for sfermion masses in the vicinity of 10^ GeV, the order of 
supersymmetry breaking scale amount required from the gauge hierarchy solution. Now, given the 
analog structure of the Yukawa interactions ([T]) and the trilinear RPV terms ^ , one can assume that 
the same (flavor-)structure responsible of the y^'"^ suppression (w.r.t. y* for the top quark) would 
also characterize the X^'jj^ couplings, inducing in turn the wanted tiny neutrino mass scale through 
one- loop processes. 

Assuming such low values for the X^jp. coupling constants (which in general easily pass the various 
phenomenological constraints [40] ) . one finds from formulas obtained in [41j that a gravitino (G) LSP 
would decay into the three SM fermion final state fifjfk via the A^^^ interactions with a sufficiently 
small width. Quantitatively, for a gravitino mass around ~ 10^ GeV, its lifetime would lie in the 
range 10^^ — 10^^ sec which is well above the age of the universe: to — 3.2 10^'^ sec. 

As a first conclusion, it is interesting to note that a scenario with trilinear RPV interactions of type A^^^ 
addressing simultaneously the neutrino mass and dark matter problems is conceivable. In the following, 
we will address this double problem within a different supersymmetric scenario where the neutrino 
mass is generated via bilinear and trilinear RPV couplings of type /Uj and Aj, respectively, which do 
not induce too large decay widths of the dark matter candidate. In contrast with the A^^.^^ case, the 
present framework deserves a more general treatment in the sense that the neutrinos are mixed with 
neutral gauginos/higgsinos so that the neutrino constraints and the neutralino ones are correlated. In 
fact, the mixing is generally so small that it does not affect significantly the gaugino/higgsino mass 
matrix, but the precise values of the several parameters entering this matrix crucially determine the 
small induced neutrino masses. 

Nevertheless, let us finish this part by commenting on the fact that even if the A-^-^^ couplings of the 
order of weak Yukawa couplings cannot induce dangerous G decay channels, they could contribute 
partially to the neutrino masses through one- loop diagrams. Here, we concentrate on the dominant 
tree-level contributions to the neutrino masses and hence leave for the future a precise and complete 
fit of neutrino data at loop-level. 

II. 3 Neutralino mass matrix 

The neutralino mass terms read. 



(4) 





(5) 



5 



in the basis "^^^ 



Kd (S) and u, 



= {B^, W^, h^, h^, s, Vi), where ^ (s) are the fermionic components of the superfields 
I = 1,2,3] denote the neutrinos. In Eq. ([5]), the neutrahno mass matrix is given, in 



a generic basis (where (z^j) = / 0, /Xj 7^ and A, 7^ 0, as will be discussed later), by 

''nms 



-Mnmssm 

03x3 



(6) 



where A^nmssm is the neutralino mass matrix corresponding to the NMSSM. For the latter mass 



matrix, we assume Vi ^ Vu,d, so that v"^ = v^ + v'^ + Si=i ''^i = 
matrix below, s and c stand for sine and cosine, respectively. 



(175GeV)2. In the mass 



-Mnmssm 



V 



Ml 


-mz Sg^ 0/3 
mz Sg^^, 3/3 







M2 

mz Cg^ Cp 
-mz Cg^ Sp 




-mz Sg^ Cj3 
mz Cg^ Cj3 



-/^ 

-\Vu 



mz Sg^ 3/3 
-mz Cg^ 3j3 








—Xv, 



2(k(s)+/x') J 



Above, Ml (M2) is the soft supersymmetry breaking mass of the bino (wino), tan/3 = Vu/vd = 
{hu)/{^d) (■^ith C/3 = COS /3 and 3/3 = sin/3), and fi = A(s). The fx' term appears only within the 
GMSB framework [see Eq.©]. 

We assume for simplicity that A, k and the soft supersymmetry breaking parameters are all real. 
In Eq. ([6]), is the RPV part of the matrix mixing neutrinos and neutralinos: 

^ ^ Q ^i + Ai(^) Ait;„ \ 



y/2 ^/2 
g'v'2 gv2 

g'vz gvs 



(8) 



V-^ ^ ^3 + A3(.) A3.„/ 



g and g' being the SU(2) and U(l) gauge couplings. We restrict ourselves to the situation where 
'^i/vu4 ^ 1 (as before), ^ 1 and |Aj/A| ^ 1 so that (i) no considerable modifications of the 

NMSSM scalar potential are induced by the additional bilinear and trilinear terms in superpotential 
dS]), (ii) the neutralino- neutrino mixing is suppressed, leading to sufficiently small neutrino masses; it 
is remarkable, as mentioned above, that the necessary order of magnitude, Vi/vd ~ Mj/m ~ Aj/A ~ 
10~^ — 10~^, corresponds typically to the hierarchy between the electron mass and the top quark mass 
(or equivalently the EW symmetry breaking scale). 



II. 4 The various (s)lepton mixings 

Chargino/neutralino-lepton mixing: The condition for generating two non- vanishing and non- 
degenerate neutrino mass eigenvalues at tree level is to ensure simultaneously 

Hi / and Ai / 0, (9) 

where Aj = (s)(Ai + A^). The respective roles of these effective parameters is reflected e.g. in the 
sub-matrix ([8]) mixing the neutrinos with neutral NMSSM states. 



6 



Since the and Lj superfields possess the same gauge quantum numbers, one can freely rotate 
La = {Hd, Li) [a = 0, ... ,3] through 5C/(4) matrices by a redefinition of fields (for examples of uni- 
tary matrix associated to SU{4:) transformations, see Ref. [l2] or Ref. [38] in Section 2.1.4). Motivated 
by the condition ([9]), we work in a general basis where the ^iLiH^ terms are non- vanishing. One could 
choose a specific basis of L^ superfields where the Aj couplings vanish, but then in general (i.e. as- 
suming no particular correlation between the Aj and Bi, rh'^- values) some Vi are generated due to the 
destabilization of the scalar potential by terms linear in originating from the RPV soft scalar bilinear 
terms ih\Li H after translation of the Higgs fields: ^ + v^j \/2 (/i^ — )• -|- v^j \Fi)- 

Reciprocally, a field basis where the Vi vanish generally leads to non-zero Aj couplings. In conclusion, 
the condition ([9]) is generally fulfilled. 

Higgs-slepton mixing: The physical amount of Higgs-slepton mixing is parametrized by the basis- 
independent angle between the 4-vectors = {vd,Vi) and Ba = {Bii,Bi), where B^ is the soft 
parameter entering the biscalar term B^huhd [38]. We assume this angle sufficiently small so that the 
additional potentially dangerous decay channels G, Xi H^t^ / hPu and u — >■ //, VV [f = fermion, 
V = Vector boson] play no role in the present context of a long-lived LSP dark matter candidate of 
type G, Xi or 9. 

Charged lepton mixing: The charged leptons also mix with the Wino (via the Wj's) and the charged 
higgsino (controlled similarly by the /ij and Aj parameters). However the detail of this mixing also 
relies on the determination of the precise value of each Yukawa coupling constant Y^ - for the charged 
leptons [see Eq.([T])] which requires a complete scenario of fiavor taking into account all the experi- 
mental constraints [from neutrino oscillations and lepton fiavor violating reactions] on leptonic mixing 
angles {Upmns matrix [l3]), a task beyond the scope of the present work. Therefore, we will not 
explicitly work out the chargino-charged lepton mixing which is expected to be of comparable amount 
as the neutralino- neutrino mixing involving similar matrix structures and parameter orders j24j . Con- 
cerning the present cosmological context, this implies for instance that the width for the decay channel 
G W^i^ is close to the one for G — t- Z^i^, so that our principal conclusions on the relative stability 
of the gravitino LSP based on the latter channel will not be modified by the former one El- 

III Sneutrino LSP 
III.l RPV sneutrino decays 

In case of a sneutrino LSP decaying through the RPV couplings of type Aj and fii, the sneutrino can 
exclusively decay into neutrinos via the two types of Feynman diagrams drawn in Fig.([T]) or into two 
charged leptons. The associated decay widths are given below. 

For that purpose, we need some definitions. First, the sneutrino LSP, noted z>i, corresponds to the 
lightest of the three sneutrino mass eigenstates, keeping in mind that arbitrary non-universal conditions 
on scalar soft masses and the dependence of sneutrino mass running e.g. on flavor-dependent Aj 
parameters lead generally to a non-degeneracy of rric^^ {i = 1,2,3) eigenvalues. 

^Similarly, there exists a trilinear RPV soft term, arising in the NMSSM, which involves the singlet scalar field; this 
one is of the form sh^Li in the Lagrangian [2^. Its direct efi'ect on the present study will be discussed in Section rill. 21 

■^Besides, note that the contributions to the charged lepton masses should be of the order of magnitude of the neutrino 
mass scale and in turn negligible compared to the direct Yukawa contributions. 



7 



We also need to introduce the matrix N^^p which is defined as fohows. For convenience, we change 
the order (w.r.t. Section III.3P of the fields in the weak basis [for the rest of the paper]: = 
{B^ , 14^3 , /i^, /i^, s, z/j) where as before the index z = 1, 2, 3, used for compact notations, corresponds to 
the three flavor states: z/g, i^^ and v^- respectively. After diagonalization, the eigenstates are ordered in 
the mass basis according to 



where v'^ denotes the three neutrino mass 



eigenstates (j = 1, 2, 3) and the five x s are the NMSSM neutralinos. Now, the unitary transformation 
matrix Nap acts as = Nap'^^ with a, /? = {1, . . . , 8}. 





Figure 1: Feynman diagrams for the sneutrino decay channel in two neutrinos, v ^ ^ 
the effective quantity rrii = fii + Xi{s) to parametrize a direct /ij mixing term effect combined with a 
(s) vev insertion. Vu denotes the up Higgs vev insertion and Vi symbolizes the sneutrino vev insertion 
(i = 1,2,3). 



Sneutrino decay into two neutrinos v^v^ (phase space factors involving neutrino over sneutrino 
masses are neglected): 



r(f>i 



m 



IQtt 



k=l 



ifc 



^/2cos9w 



(10) 



where i,j = {1,2,3} run over the neutrino eigenstate indexes, Mik are elements of the model- 
dependent sneutrino basis transformation matrix (1 is for the lightest sneutrino eigenstate and 
k = {1,2,3} is for the sum over sneutrino flavor states), U = —NnsinOw + Ni2Cos6w and 

i 

o j]' indicates that the same expression in brackets must be added by switching i with j. 

Sneutrino decay into two charged leptons C-fi^ (phase space factors involving lepton over sneu- 
trino masses are also negligible here): 



m, 



IQtt 



k=l 



Ik 



+ 5.. 



Uj{k+2}Ue^w 



in) 



where C/j|fc_|_2} is the equivalent of neutralino rotation matrix but for charginos (here with 
k = {1,2,3} for the sum over charged lepton flavor states and j = {1,2,3} associated to the 
final charged lepton eigenstate) and similarly U^,^^ is the equivalent for charginos of U^m^- 
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III. 2 Sneutrino stability 

It turns out that the strength of effective RPV couphngs Aj/(s) and /Xj necessary to generate sufficiently 
large neutrino masses, i.e. typically m^-^ > 10^^ — 10^^ eV (in the approximate case of a unique 
neutrino flavor), induces RPV-like mixings between neutrinos and neutralinos which are too large 
from the cosmological point of view. Indeed, those mixings translate into matrix elements Ni^ ~ 10"'', 
Nj5 ~ 10"''', Uj^m^ ~ 10"'^ and Nj^j^_^^j ~ 1 with i,j,k = 1,2,3 [for standard neutralino masses 

m^o = 0(10^) GeV] so that the total sneutrino RPV decay width T{i>i z^f ) ~ 10"^^ GeV is well 
above the critical value of 1.52 10~^^ GeV [for a sneutrino mass: > 10 GeV and assuming for 
now A/ifc ~ 1]. It means that the sneutrino lifetime is well below the age of the universe and thus a 
sneutrino LSP does not constitute a viable dark matter candidate. A source of width suppression may 
arise from the generic amount of the matrix element Mu but its effect is at much of a few orders of 
magnitude (more suppression would correspond either to a fine-tuning of parameters or to the weak 
breaking of a certain symmetry which should be described) so that this cannot invalidate the above 
conclusion. Furthermore, the effect of the additional decay channels into two charged leptons and into 
(pseudo) scalar Higgs bosons, h^h^, h^a^ or a^a^ [i,j = 1,2,3, k,l = 2,3] (via the soft term shuLi 
which is characteristic of the NMSSM with RPV interactions), can only increase the total sneutrino 
width 0. 

IV Neutralino LSP 

IV. 1 RPV neutralino decays 

If the LSP is the lightest neutralino and the RPV couplings of type A, and are present, then the 
lightest neutralino can only decay into the W,Z bosons or into the (pseudo) scalar Higgs fields [as 
illustrated by the Feynman diagrams of Fig. ([2])- ([3])] if one restricts oneself to the dominant two-body 
decay channels. The obtained partial decay widths are given as follows. 




Figure 2: Feynman diagram for the neutralino decay mode into the Z boson, Xi — ^ Z uj^. The cross 
allows to specify the neutralino component that is coupled (here: a neutrino flavor state Vi). 



''in this part, we study the sneutrino RPV decays for completeness, even if a left-handed sneutrino LSP as a candidate 
for dark matter has been already excluded by direct dark matter searches for most of the realistic sneutrino mass ranges. 
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Figure 3: Feynman diagram for the neutralino decay modes into scalar Higgs fields, Xi ~^ ^k'^j^^ 
pseudoscalar Higgs fields, x? — ^ The cross indicates which neutralino components are coupled. 



Neutralino decay into a Z boson and a neutrino v^' (phase space factors involving neutrino 
masses are neglected): 



1 

96^ 



i=l 



T7^^y° 



1 



1 + 



2m| 



(12) 



where j = 1,2,3 labels the neutrino mass eigenstate and z = 6, 7, 8 corresponds to the sum over 
neutrino flavor states. 

Neutralino decay into the W boson and charged lepton (the dependency on lepton masses, 
which represents subleading effects, is omitted): 



r(x? 



1 

48^ 



i=6 



m 



w 



2 



1 + 



Xi 



(13) 



where j = 1,2,3 labels the charged lepton mass eigenstate and i = 6, 7, 8 is corresponding to the 
charged lepton/neutrino flavor states. 



Neutralino decay into scalar Higgs bosons /i? and neutrino z^^: 



r(x? ^ hluj^) = — 

ozvr 



8 
i=6 



|A'435'a:3 + ^455'fei|^m,-o 1 



2 



(14) 



with J = 1,2,3 labeling the neutrino mass eigenstate, i = 6,7, 8 corresponding to the neutrino 
flavor states and /i^ = /i^i ^3 (scalar Higgs mass eigenstates). The rotation matrix 5 relates 
the real parts of the neutral Higgs bosons and singlet scalar field to the scalar Higgs mass 
eigenstates /i^ (see the precise definition in ^^j). We recall that the parameters A|j_5} which 
appear here are the NMSSM specific trilinear parameters of Eq.©. We also mention the decays 
into charged Higgs fields which are expected to be of comparable widths. 

Neutralino decay into pseudoscalar Higgs boson and neutrino v"^: 



r(x? 



1 



{i-5} 



i=6 



|iV43^fc3 + A^45AlP"l-o 1 




(15) 
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with j = 1,2,3 labeling the neutrino mass eigenstate, i = 6,7,8 corresponding to the neutrino 
flavor states and = a?, 02 (pseudoscalar Higgs mass eigenstates). The rotation matrix P 
translates the imaginary parts of the neutral Higgs bosons and singlet scalar field into the 
pseudoscalar Higgs mass eigenstates [16]. 

IV. 2 Neutralino stability 

Based on the typical values of RPV parameters A, , Vi and //j generating a neutrino mass scale niui ~ 
10~^ — 10~^ eV, we find for the parameters involved in the partial widths (jl2p . (jl3p . (|14p and <\15\i : 
Nj, ~ 1, ~ 10-^ [/,|,_3| ~ 1, A|,_5| ~ 10-6, iV435fc3 + Ni^Ski ~ 10"! - lO'^ and A^43Pfc3 + 
~ 10-^ - 10-2 with i = 6,7,8 and j = 1,2,3 [if neutralino masses are of the order: m^o = 
0(10^) GeV, and taking m^o > mz]- We can see that the contributions to the total neutralino 
decay width of Eqs. (jl4p and (jlSp are smaller than Eqs.(|12p and (jl3p . This leads [whatever is the 
(pseudo) scalar Higgs spectrum] to a total neutralino RPV decay width r^otai (X?) ~ 10-15 GeV that 
is several orders of magnitude above the critical value of 1.52 10"^^ GeV. There are even additional 
contributions e.g. to the lightest neutralino decay into Z^uJ^, of the same order (originating from the 
Z^XiX^ coupling), that should slightly increase the total neutralino decay width. Adding the decay 
channel into H^, closed in most of the parts of the parameter space, can only increase again the total 
neutralino width. The conclusion is thus as for the sneutrino case: the neutralino lifetime is much 
smaller than Iq and hence it does not represent a possible dark matter LSP candidate. 

Nevertheless, a possibility to insure the stability of the lightest neutralino is to restrict to domains 
of parameter space where the 4 possible decay channels of Eq.([12|), Eq. (fT3]) . Eq. p^ and Eq. p^ . as 
well as the channels into charged Higgs bosons, are all kinematically closed [even opening only one 
channel is sufficient to render the Xi unstable]. Although in the NMSSM the lightest pseudoscar Higgs 
boson noted ai can be much lighter than in the MSSM, we can find realistic regions of the NMSSM 
parameter space where the lightest neutralino can be simultaneously lighter than the W, Z bosons, 
the lightest scalar and pseudoscalar a^, as well as the charged Higgs bosons H^, thus closing 
the corresponding RPV channels. By the term 'realistic', we mean that the NMSSM parameters 
pass the theoretical and phenomenological constraints implemented in the NMHDECAY program [16] 
like: (i) the physical minimum of the scalar potential is deeper than the local unphysical minima 
with {h^^) = and/or (s) = (ii) the running couplings A, k, Y^, do not encounter a Landau 
singularity (in) the experimental constraints from LEP in the neutralino, chargino and Higgs sectors 
are effectively satisfied. 

For instance, using NMHDECAY, we have plotted the masses m^o , m^jO and m„o as a function of the 
K coupling constant in Fig. ([4]), for fixed values of the other NMSSM parameters. We see on this figure 
that for K > 0.05, the lightest neutralino LSP is lighter than the W, Z bosons, the lightest scalar 
and pseudoscalar a^. So the scenario of a completely stable Xi LSP can be a priori realized. Note that 
in such a scenario the RPV couplings can now be chosen freely to satisfy the neutrino constraints, 
since those couplings have a negligible impact on the neutralino mass spectrum. 

We remark that there exist various triangular one- loop processes, exchanging charged leptons/sleptons, 
quarks /squarks or charged leptons /W bosons (through bilinear or trilinear RPV couplings El) , which 
contribute all to the decay channel Xi ~^ l^T'- Such a channel is always kinematically open in 
the NMSSM, even in the case of the specific spectrum discussed just above, and could thus render 

^The contributions from trilinear RPV interactions were estimated in Ref. [44| . 
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Figure 4: Masses m^o [purple curve], m^jO [green curve] and m„o [red curve] as functions of k. We 
fix the other NMSSm' parameters: tan/3=2, A=0.7, /i = A(s)=53d GeV, Mi=66 GeV, M2=133 GeV, 
M3=500 GeV, >1a=1280 GeV, ^^^=0, At = Af, = Ar=-2.b TeV, m^-±=200 GeV and mq=l TeV 
(universally). The A parameters are the trilinear scalar soft supersymmetry breaking parameters 
which do not affect the neutralino mass matrix. The low values of Mj's allow to get a neutralino LSP. 
Note also the particularly low values of m^o, characteristic of the NMSSM. For completeness, we give 
the employed RPV parameters: Ai/(s) = 9 10~^ and //i = 6 10~^ GeV. 



the neutralino unstable from the cosmological scale point of view. Assuming that no important 
destructive interferences occur among the several triangular contributions mentioned above, rough 
estimates indicate that the suppression provided by RPV terms [reproducing the correct neutrino 
mass scale] and by loop factors are not sufficient to make the LSP neutralino stable with respect to 
to (whatever is the neutralino composition). However, precise loop calculations of all contributing 
reactions, beyond the scope of our study, are necessary to conclude definitively on this aspect, taking 
into account notably effects of heavy sfermion loop factors. 

The other important final comment is that three-body decay channels like Xi ~^ £'^£~u"^ (via an 
off-shell Z/W boson [45^ 146] or an off-shell Higgs boson 05]) are systematically open. From order of 
magnitude estimates, it tuns out that such channels are expected to render the lightest neutralino 
clearly unstable. 
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V Gravitino LSP 



V.l RPV gravitino decays 

In the case of a gravitino LSP and RPV couplings of type Aj and /ij, the gravitino can decay into the 
EW gauge bosons and the (pseudo)scalar Higgs fields [as shown in Fig.(l5])-([6])]. We give below the 
obtained associated partial decay widths. 

All formulas given for the gravitino width, and above for sneutrino and neutralino widths, result 
from original calculations. Nevertheless, some of the gravitino decay amplitudes, into 7, Z^, 
and , were computed with other conventions by the authors of Ref. \20\ [23| ITT] in the case of the 
MSSM (see also Ref. for a more generic calculation approach). Note in particular that within the 
present framework of the NMSSM, the numerical results for the gravitino decay into the pseudoscalar 
- which can be quite light there - might differ significantly from the MSSM and thus decrease 
dangerously the gravitino lifetime. Moreover, within the NMSSM, new contributions arise for the 
decays G — )■ /i^ff^ and G — a^i'l^ due to the presence of the s superfield [see the right Feynman 
diagram of Fig. ([6])]. Finally, in the NMSSM, the neutralino mass matrix (and in turn the matrix 
elements Nij involved in gravitino amplitudes) takes a specific form, especially in the context of the 
GMSB (see Section HlHI) . 




Figure 5: Feynman diagrams for the gravitino decay processes into a photon, G — ?■ 7f[", and a Z 
boson, G —7- Z^u^. Vi denotes a sneutrino vev insertion. 



Gravitino decay into the photon 7 and neutrino - the dependency on the neutrino mass is 
omitted (see first calculation in Ref. |5]): 

where Uu^^'^ = NucosOw + Ni2sin9w, i = 1,2,3 labels the neutrino mass eigenstate and 1 (2) 
corresponds to the (W^) component. Mpi is the reduced Planck mass: Mpi ~ 2.4 10^^ GeV. 

Gravitino decay into the Z boson and neutrino (see original calculation in Ref. 

r(6^z»o = J^|j/„,,|^;|j(i-^| 11 + .^ + .^ I ,17) 
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Figure 6: Feynman diagrams for the gravitino decay processes into scalar Higgs fields, G — t- h^vf^^ 
and pseudoscalar Higgs fields, G — )• a^vf^- We use the symbolic notation rrij = /Xj + Aj(s) combining 
a direct /Xj mixing and a {s) insertion. Vu indicates an up Higgs vev insertion. 



where i = 1,2,3 labels the neutrino mass eigenstate and 1 (2) corresponds to the (W^s) 
component. 

Gravitino decay into the W boson and charged lepton tf (see original calculation in Ref. |24j): 

no WH-) = ^\u,^\'^ f 1 - f 1 + ^4^ + (18) 

where numerically U^,-^^ is expected to be of a comparable order to U^m^ as already discussed. 
Similarly, here, i = 1,2,3 labels the charged lepton mass eigenstate. 

Gravitino decay into the scalar Higgs boson and neutrino 

1 'ml/0 ( '^lo \^ 

3847r y mli^) 

with i = 1, 2, 3 labeling the neutrino mass eigenstate, the numbers 3, 4, 5 corresponding respec- 
tively to the h^, hj^, s components and 1, 2, 3 respectively to the real h^, h^, s components. 

Gravitino decay into the pseudoscalar Higgs boson and neutrino uf^ 

still with k= 1,2. Similarly, i = 1,2,3 labels the neutrino mass eigenstate, the numbers 3, 4, 
5 correspond respectively to the /i^, /ij], s components and 1, 2, 3 respectively to the imaginary 
/i^, /ij], s components. 

Gravitino decay into the charged Higgs boson and charged lepton tf 
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i = 1, 2, 3 labels the charged lepton mass eigenstate. Ui2 represents the components of the three 
charged lepton eigenstates if into the charged higgsino h^. 

We note that the partial widths for the charge conjugated final states are equal, which means for 
instance that T{G — t- ji)) = T{G — t- 71^). Later, we will thus refer to the total gravitino decay 
width as: rtotai(G) = 2[r(G ^ 71^™) + T{G Z^u"") + T{G W+f-) + T{G h^u"') + T{G 
a^u'^) + T{G ^ H+r)]. 

V.2 Gravitino stability 

V.2.1 Scenario I: almost all decay channels open 

Let us first consider a NMSSM scenario which is a priori dangerous from the cosmological point of 
view: we choose an heavy gravitino, so that the 5 first of the 6 types of decay channels described 
above are kinematically open (parameters allowing only the decays into the first h\ and o^), and a 
quite light pseudoscalar a\ (as is possible in the NMSSM), tending to increase the phase space for the 
partial width T{G a^z^f*). 

First, the considered region of the parameter space must be such that the gravitino is the LSP. We 
thus take rather large values for the gaugino masses Mi, M2 and M3: Mi=300 GeV, M2=600 GeV, 
M3=2 TeV to push the neutralino mass to higher values. For the same reason, we take m^±=300 GeV 
and mq=l TeV. We fix the other NMSSM parameters at: A=0.3, /i=237 GeV, «;=0.35, A«:=-30 GeV, 
Af = A), = j4t-=-2500 GeV (trilinear soft parameters). We allow tan/3 and Mj^ to vary accordingly to 
a scan performed for tan/3 = 4—7-14, Ma = 250 — t- 300 GeV. Ma represents the pseudoscalar mass 
in the MSSM but it is an effective parameter in the NMSSM. It is somewhat equivalent to the second 
(pseudo) scalar masses in the NMSSM and to the charged Higgs mass which are almost degenerated in 
this model: Ma — mi^o,m^o,mjj±. It is related to the other ones via the minimization equations. In 
the Fortran code NMHDECAY, we can choose between Ma or A\ as input parameters. This scan is 
performed through the NMHDECAY code so that the generic NMSSM constraints mentioned above 
are satisfied for the selected parameters. 

As a first step, we present in Fig.© the output masses m^o, mc,^, m^o and m^o, obtained via this 
scan, as functions of tan/3. 

Among the possible points of Fig.©, we choose the one corresponding to tan/3=8 and M^=275 GeV 
leading to m^o= 219.5 GeV, mi>=293.5 GeV, ^1^0= 116.7 GeV, m^o= 259.9 GeV, mj,Q= 548.8 GeV, 
m„o= 25.3 GeV, 1^0°— 303.5 GeV and mfj±= 265.5 GeV. The gravitino mass is fixed at 7713/2=200 GeV 

so that G is well the LSP. Then the channels G hlv'^, G alv"^ [with k>2] and T{G H^t^) 
are kinematically closed but all the other ones are open. 

Now, in the one lepton flavor approximation we choose the RPV parameter values Ai/(s) = 2 10~^ 
and /Ui = 10~^ GeV leading to m^^ = 1.04 10^^^ GeV^ which is reasonable from the point of view of 
experimental neutrino data. 

Finally, for the chosen NMSSM parameters and RPV couplings, the induced RPV neutrino-neutralino 
mixings give rise to the partial widths T{G 71^") = 1.95 lO"''^ GeV, T{G Z'^u"^) = 8.80 10"^^ 
GeV, r(G ^ W+r) = 1.94 10"^^ GeV, T{G /i?z/"^) = 9.47 10"^^ GeV and r{G a?i/™) = 
2.08 10~^^ GeV. It is important to note that the relative smallness of T{G — )• 72^*") reflects in particular 
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Figure 7: Masses m^o [purple points], mj> (universal) [blue points], m/jO [red points] and m^o [green 
points] as functions of tan^. We use: A=0.3, k=0.35, //=237 GeV, Mi=300 GeV, M2=600 GeV, M3=2 
TeV, ^«;=-30 GeV, At = Ah = ^^=-2500 GeV, m^-±=300 GeV (universal), mq=l TeV. The points are 
obtained from a scan performed on the two parameters tan/3 = 4 — )• 14 and Ma = 250 — t- 300 GeV, 
using the NMHDECAY code [H]. The chosen RPV parameters are: Ai/(s) = 2 10"^ and m = 10"^ 
GeV. 



the smallness of the photino component for the neutrino mass eigenstate. The corresponding total 
gravitino width is FtotaKG) = 2.00 10~^^ GeV giving rise to a gravitino lifetime Tq ~ 76.3 Iq. This 
result illustrates the feature that within the NMSSM, and for RPV couplings that generate the correct 
neutrino mass scale, a gravitino LSP is stable (with respect to to) even for a gravitino mass up to 200 
GeV typically in agreement with supergravity (or mixed with GMSB) scenarios. 
Indeed, more generally speaking F^otal 

(G) remains of the same order of magnitude as for the above 
NMSSM parameters, if rrii, ~ 10""^"'^ GeV (as imposed by data) and m^o ~ '"1.3/2 ~ 10^ GeV (as 
natural in supergravity-like models). The reason, taking e.g. F(G — t- 71^"^), is that the orders of rrii, 
and m-j^o fix the order of Ui,m;y (see Eq.(6) in Ref. [5j) and in turn the partial gravitino decay width 
(c./. Eq. (fT6]) ) once m3/2 is chosen. Furthermore one expects C/^m^ — U^^^ ~ f/j^m^ in Eq. (fT7|) - (fT8]) 
[21]. Similarly, the partial decay widths into the (lightest) neutral Higgs bosons are fixed by m^o, m^o, 
m^l2 as well as the amount of neutrino components in the higgsino and singlino (c./. Fig.©) [and 
thus by my and itl^o\. If kinematically open, the decays into the h2/a^ Higgs bosons have rates of 
similar order (for identical parameters). For a a\ boson heavier than in the above set of parameters, 
T{G ^ a^u"^) decreases but the total width remains around ~ 10 GeV. Finally, if allowed, the 
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channel G — )• H^l^ is expected to reach a width of the same order as G — )• h^v^ and hence would 
not significantly modify Tiota\{G)- 

It has been recently 09] pointed out that the widths for G (.^W^"- i^ff and G u'^Z* 
i/^ff together can reach 10^ x T{G — )• 71/™) restricting oneself around parameter domains where 
m^o ~ m3/2 ~ 10^ GeVlfl. It means that the new three-body decay widths can reach ~ 10~^^ GeV, 
from our present results, and this does not modify rtotal(G) significantly. 

V.2.2 Scenario II: a few decay channels open 

To be more general than in the approach of Section rV.2.1l we consider here the case where the gravitino 
is the LSP and is lighter than the W boson as well as the (pseudo)scalar Higgs fields [rather than 
considering some given points of the parameter space]. Then the RPV channel through the process 
G — )• 71^™ remains open. Imposing the width r(G — ?■ 71^™) to be smaller than the critical value of 
1.52 10~^^ GeV leads to the upper bound on the gravitino mass ^ 3 TeV for \U^u\ ~ 10~^ (see 
Eq.([T6|)) as imposed by the neutrino mass scale if m^o < 10^ GeV. For example, let us consider a 
given point of parameter space for which such a neutralino mass is realized. We take randomly the 
point of Fig.© with k=0.1. Then, m^o= 67.8 GeV, m^=193.9 GeV (universal), m^o= 128.7 GeV 
and m^o= 83.2 GeV [insuring also that the gravitino decays into the (pseudo)scalar Higgs fields are 
forbidden]. For this point, one gets Uuy^^ = 1.94 lO"'^ with Ai/(s) = 9 10"'^ and m = 6 10"^ GeV 
which reproduce the correct neutrino mass scale (m^^ = 5.73 10~^^ GeV^). The exact bound resulting 
from the gravitino lifetime in this precise case is 771,3/2 < 2.7 TeV. Taking into account the other allowed 
channels, G — )■ i^ff' and G — > v^ff [/ = fermions] (see Ref. [49]), this bound changes to 777.3/2 ^ 650 
GeV - 3 TeV as long as 777^0 ~ 10^ GeV. Since we are in a situation where 7773/2 < mw, this bound 
is respected. From the general point of view, we conclude that the NMSSM with phenomenologically 
interesting RPV couplings implies a gravitino LSP which is systematically stable (with respect to to) 
if one takes its mass to be smaller than myi/, m^o and 777^0. 

This quite general conclusion applies to supergravity (or mixed with GMSB) scenarios where 7773/2 ~ 50 
GeV is realistic, as well as to pure GMSB scenarios in which possibly 7773/2 ~ 1 eV. Indeed, within the 
GMSB framework, the new form of the neutralino mass matrix (see Section III.3P leads to values of 
Unv^Xf extremely close to the above case [and hence to same bounds on 7773/2] for the realistic neutrino 
mass scale, if one chooses e.g. the typical characteristic value ^' = 100 GeV. Note that in GMSB, the 
gravitino is clearly automatically lighter than W , and a^. Note also that the channel G — t- 77^^" is 
systematically open (in the realistic 3 flavor case) since within our neutrino scenario rriu^ = at the 
tree level. 

V.2.3 Scenario II extended to 3 neutrino flavors 

There exist sets of NMSSM parameters [passing the NMHDECAY constraints] and RPV couplings 
which reproduce the squared neutrino mass eigenvalue differences measured in oscillation experiments. 
More precisely, the updated three-flavor analyzes based on a global fit including results from solar, 
atmospheric, reactor and accelerator oscillation experiments lead to (at the 3(T level): 7.1 < Am^]^ < 
8.3 [lO-^eV^] and 2.0 < A777|^ < 2.8 [lO'^eV^] [13] . In our notations e.g. A777I1 = 1777^^ -777^ J. Values 

^This ratio result is not expected to be significantly changed when extending the MSSM to the NMSSM. 
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lying in those two intervals arise for instance for the point of Fig.dl]) if k=0.1 and0 

^1 = 3 10^^ GeV, fi2 = 1.5 10-^ GeV, ^3 = 2.5 10"^ GeV 

Ai/(s) = 8.5 10-^ A2/(s) = 1 10-^ A3/(s) = 1.5 10"^. (22) 

Indeed, this complete set of parameters yields the following three neutrino mass eigenvalues at tree 
level: 

m^^ = eV, = 0.00867 eV, m^^ = 0.04670 eV (23) 

in agreement with the above experimental intervals for Am2i and Am-|^ since in the present RPV 
model one has Amgi = m^^ = '''•^l 10^^ eV^ and Am^i = = 2.18 lO^'^ eV^. 
Prom the cosmological point of view, the neutrino mass eigenvalues in (I23p satisfy the bound extracted 
from WMAP and 2dFGRS galaxy survey (depending on cosmological priors): X]i=i "^^-i ~ ^■'^ 

m- 

Finally, the above neutrino mass eigenvalues are perfectly compatible with the limits extracted from 
the tritium beta decay experiments (95% C.L.): mp < 2.2 eV [Mainz] and < 2.5 eV [Troitsk], this 
effective mass being defined as = X^^^i \ Uei\'^rn'^. where Uei is the leptonic mixing matrix [5JJ. 

The point of parameter space considered in this section is the same as the one considered in the previous 
section (i.e. Section [V.2.2p . namely: the point of Fig.(|l]) with k=0.1, and we still assume a situation 
where the gravitino LSP has a mass smaller than mw, m^o and m^o (scenario II). For this point, 
and with the RPV couplings of Eq. (j22p . we obtain for the gravitino decay channel into a photon: 
C/^m^ = -4.74 10~^\ U,y^^ = -6.61 10^*^ and C/^.^™^ = 2.30 lO"'^, recalling the definition U,y^^ = 
NiicosOw + Ni2sin6w where i = 1,2,3 labels the neutrino mass eigenstate and 1 (2) corresponds to 
its (^3 ) component. Those numbers imply the limit 711-3/2 < 2.5 TeV resulting from the condition 

2 ^(G^li^r) < 1-52 10"^^ GeV. 

i 

The resulting bound < 2-5 TeV obtained here is close to the bound ^71,3/2 < 2.7 TeV of previous 
section, for the channel G — )• 7/^*", which means that moving to the three flavor case should not affect 
significantly the conclusions. Therefore we conclude that, in the case of three neutrino flavors as well, 
a gravitino LSP is always sufficiently stable (with respect to to) if its mass is weaker than mw, m^o 

and m„o. 

"1 

This conclusion on the 3 neutrino flavor case is not trivial in the sense that, starting from the 1 
flavor situation, the variation of the f/i/^^ value when extending to 3 neutrino flavors cannot be easily 
predicted due to the rich structure of the whole RPV neutralino mass matrix. Moreover, there is no 
simple argument to deduce the values of 11^^;^ and Ui,j^^ from the [7,^™^ element (encoding the neutrino- 
gaugino mixing) because of the multiple mixing types between the different ffavors of neutrinos and 
the various neutralino states. 

Nevertheless, we have checked that, as expected, one obtains numerically the following hierarchy 
among the following matrix elements: |A'^ij| < |-^2j| < |-^3j| [j = 1)2]. This is interpreted by the fact 
that here the neutrino eigenstate is heavier than (and in turn rriu^ > rrii,^) since the massive 
neutral gaugino [B^ , VF3] components of v'!^ are larger than in the z^™ state (in turn, than in z^{"). 

^Let us note that we have chosen the specific set ((22} of RPV parameters to iiiustrate on an expiicit exampie that 
the three fiavor neutrino data can be reproduced. Nevertheless, the orders of magnitude in Ea. (|22p of the six effective 
RPV parameters are general in the sense that those are imposed systematically by the experimental ranges for Am2i 
and Amli as soon as ra^ o ^ 10^ GeV (as occurs in supergravity and GMSB models). Hence the conclusions at the 
end of this Section IV. 2. 31 which are based on these orders of magnitude, can also be generalized to any region of the 
parameter space where m^o ~ lO'^ GeV. 
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VI Discussion and conclusions 



In the context of the NMSSM with the presence of RPV couphngs large enough to generate reahstic 
neutrino masses, a gravitino LSP of mass 0(10^) GeV ~ as appears in supergravity models ~ is 
sufficiently stable from the point of view of the age of the universe. 

Nevertheless, the gravitino lifetime that we obtain is of order ~ 10^^ sec. Now even in a supergravity 
situation with 771,3/2 ^ 80 GeV, where the opened decay channels G — )• 71^™, G — )■ ff and G — )• ff 
are thus reduced, the total width looses typically four orders of magnitude only. Hence, the HEAT 
excess in the positron fraction [7J and the exotic positron source apparently detected by PAMELA 
[20] (which both require, in standard solutions, ^3/2 ~ 100 — 200 GeV and Tq ~ 10^^ sec) do not seem 
to be simultaneously explainable by the present RPV scenario reproducing the neutrino masses. The 
reason is that the neutrino-neutralino mixing necessary to create large enough neutrino masses seems 
to induce too large gravitino RPV decay widths. 

Moreover, this result that the gravitino lifetime is systematically smaller than the lifetime needed 
to explain the PAMELA excess with a gravitino dark matter also means that, in our dark matter 
scenario, the gravitino decays produce always too large fluxes which are excluded by PAMELA in 
particular. Hence our conclusion is that a gravitino LSP cannot be a good dark matter candidate if 
there exist significant RPV mixing terms (reproducing the neutrino masses). 

In the case of the existence of a special supergravity scenario ~ like e.g. some hybrid supergravity- 
GMSB models ~ where one could have 013/2 below ~ 10 GeV (in order to sufficiently reduce r(G — t- 
71/™)) and Ml ~ M2 ~ 10^ GeV (suppressing the T{G i^ff',v'^ff) contribution f39]), the 
gravitino lifetime would then be above the PAMELA limit of ~ 10^^ sec, making G a stable and viable 
dark matter candidate [with RPV couplings producing the neutrino masses]. 

Now if the supersymmetry breaking relies instead on a pure GMSB mechanism, the neutralino mass 
matrix includes new elements and the gravitino mass can decrease typically to eV scale. Then, the 
gravitino is clearly a stable LSP since we obtain a bound 7713/2 ^ 650 GeV - 3 TeV from the requirement 
that its lifetime exceeds the age of the universe, if the photino component of the neutrino [of value 
comparable to the above supergravity case] is sufficiently large to induce correct neutrino masses. 
In this GMSB case, the gravitino lifetime is ~ 10^^ sec. Hence the PAMELA data cannot be explained 
with a gravitino dark matter. However, the PAMELA flux excess can have e.g. alternative natural 
and astrophysical explanations, like in electron-positron pairs produced by nearby pulsars [52] By 
consequence, the GMSB framework could allow for the existence of a good gravitino dark matter 
candidate within the RPV-NMSSM. 

Finally, we have checked numerically the robustness of the above conclusions on the gravitino stability 
when extending to the realistic case of three neutrino flavors. To illustrate this, we have selected 
some RPV-NMSSM parameter sets in agreement with the various constraints implemented in the 
NMHDECAY code and reproducing the squared neutrino mass eigenvalue differences measured in 
oscillation experiments. 

Concerning the stability of a sneutrino LSP with respect to to) within RPV versions of the NMSSM 
reproducing the wanted neutrino mass scale, the result is negative. The lightest neutralino LSP is 

*We also note that, anyway, the mentioned interpretation [7] of the HEAT anomaly predicts an antiproton flux which 
tends to be too large, although the prediction suffers from significant uncertainties. Besides, electron and positron 
fluxes from gravitino decays (whatever are the gravitino characteristics) cannot explain |20| both the PAMELA positron 
fraction and the electron plus positron flux measured by Fermi LAT |53| . 
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also not expected to be stable with respect to to- 

Finally, we comment that all these difficulties one faces for having a viable dark matter candidate in 
RPV models are to be bridged with the astrophysical problems one encounters in baryogenesis when 
RPV interactions are indeed turned on to generate the neutrino masses [54J- 
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